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TOWARDS A THEORY OF QUANTUM COMPUTABILITY 


STEFANO GUERRINT, SIMONE MARTINI § , AND ANDREA MASINI 11 


Abstract. We propose a definition of quantum computable functions 
as mappings between superpositions of natural numbers to probability 
distributions of natural numbers. Each function is obtained as a limit 
of an infinite computation of a quantum Turing machine. The class of 
quantum computable functions is recursively enumerable, thus opening 
the door to a quantum computability theory which may follow some of 
the classical developments. 


1. Introduction 

Despite the availability of a large corpus of result^] quantum computabil¬ 
ity still lacks a general treatment akin to classical computability theory. 
Taking as a reference model (quantum) Turing machines, one of the main 
obstacles is that while it is obvious how to understand a classical Turing 
machine (TM) as a device computing a numerical function, the same is not 
so for a quantum Turing machine (QTM). 

In a naive, but suggestive way, a QTM may be described as a classical 
TM which, at any point of its computation, evolves into several different 
classical configurations, each characterised by a complex amplitude . Such 
different configurations should be imagined as simultaneously present, “in 
superposition ”—a simultaneity formally expressed as a weighted sum ^ 
of classical configurations Q, with complex coefficients d{. Even when start¬ 
ing from a classical configuration, in a QTM, there is not a single result, but 
a superposition from which we can read off several numerical “results” with 


Date : April 10, 2015. 

* Partially supported by the Project ELICA (ref. ANR-14-CE25-0005), of the ANR 
program “Fondements du numerique (DS0705) 2014”. 

^ Partially supported by the Italian “National Group for Algebraic and Geometric 
Structures, and their Applications” (GNSAGA-INDAM). 

^ Partially written while at LIPN, Institut Galilee, Universite Parisl3, Sorbonne Paris 
Cite as visiting researcher. 

1 See, in the large literature, 

nm na Hu m ns] for fundamentals results, [3] for the 
foundations of quantum complexity, or P US] |6] 0 E 0 US EH EEl El] for more language 
oriented papers. 


1 





2 


S. GUERRINI, S. MARTINI, AND A. MASINI 


certain probabilities Moreover, QTMs never have genuinely finite compu¬ 
tations and one should therefore find a way to define when and how a result 
can be read off. 

In this paper we propose a notion of “function computable by a quantum 
Turing machine,” as a mapping between superpositions of initial classical 
configuration^ to probability distributions of natural numbers, which are 
obtained (in general) as a limit of an infinite QTM computation. 

Before reaching this point, however, we must go back to the basics, and 
look to the very notion of QTM. Because, if it is true that configurations of a 
QTM are superpositions ^ g^Q of classical configurations, quantum physics 
principles impose severe constraints on the possible evolutions of such ma¬ 
chines. First, in any superposition we must have \di\ 2 3 = 1- Sec¬ 

ond, there cannot be any finite computation—we may of course name a state 
as “final” and imagine that we read the result of the computation (whatever 
this means) when the QTM enters such a state, but we must cope with 
the fact that the computation will go on afterwards. Moreover, since any 
computation of a QTM must be reversible [2], in the sense that the operator 
describing the evolution of the QTM must be invertible, we cannot neither 
force the machine to loop on its final configuration. On the other hand, be¬ 
cause of reversibility, even the initial configuration must have a predecessor. 
Summing up, an immediate consequence of all the above considerations is 
that every state must have at least one incoming and one outgoing transi¬ 
tion and that such transitions must accord to several constraints forced by 
quantum physics principles. In particular, transitions must enter the initial 
state—since a priori it might be reached as the evolution from a preceding 
configuration — and exit the final state—allowing the machine configuration 
to evolve even after it has reached the final result. 

The reversibility physical constraints are technically expressed by the re¬ 
quirement that the evolution operator of a QTM be unitary. If we now want 
to use a QTM to compute some result, we are still faced with the problem 
of when (and how, but for the moment let’s postpone this) reading such a 
final result, given that the computation evolves also from the final state, and 
that, without further constraints, it might of course evolve in many, different 
possible ways. Bernstein and Vazirani in their seminal paper [3] (from now 
on we shall refer to this paper as “B&V”) first define (non unitary) QTMs; 
select then the “well-formed” QTMs as the unitary-operator ones; and de¬ 
fine finally “well-behaved” QTMs as those which produce a superposition in 
which all classical configurations are simultaneously (and for the first time) 
in the final state. What happens after this moment, it is not the concern of 
B&V. 

2 We cannot observe the entirety of a superposition without destroying it. But if we 
insist on observing it, then we will obtain one of the classical configurations Cj, with 
probability \di \ 2 . 

3 More precisely, the domain of our functions will be the Hilbert space of square 

summable, denumerable sequences of complex numbers, with unitary norm. 
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Our goal is to relax the requirement of simultaneous “termination”, allow¬ 
ing meaningful computations to reach superpositions in which some classical 
configurations are final (and give a “result”), and some are not. Those which 
are not final, should be allowed to continue the computation, possibly reach¬ 
ing a final state later. The “final result” will then be defined as a limit of 
this process. In order to assure that at every step of the computation the 
superposition of the final configurations is a valid approximation of the limit 
result, we must guarantee that once a final state is entered, the “result” that 
we read off is not altered by the further (necessary, by unitarity) evolution. 
To obtain this, we restrict the transition function out of a final state, without 
violating unitarity. We obtain this goal by “marking” the symbols on the 
tape—once a final state is entered, the machine remains in that final state, 
replaces the symbol a under the head with the same marked symbol a, and 
moves to the right. Dually (to preserve unitarity), when the machine is in 
an initial state, it may rollback to another initial configuration in which the 
symbol a to the left of the head is replaced by the marked symbol a and the 
head is on it; that is, looking at the corresponding forward transition, when 
in the initial state, if the machine head reads a marked a, then the machine 
remains in that initial state, replaces the symbol with the unmarked a, and 
moves to the right. The role of the extra marked symbols is restricted to 
these “final” and “initial” evolutions. In particular, there are no transitions 
out of a final state when reading a marked symbol, or which enter an initial 
state writing a marked symbol, and no transitions at all—involving marked 
symbols—entering or exiting a state that is neither final nor initial. That 
these machines correctly induce a unitary operator is the content of Theo¬ 
rem [7| In the following, marked symbols will be called the extra symbols 
and we will generalise this discussion from a final (initial) state to a set of 
target ( source ) states. 

After the definition of QTMs and of the corresponding functions, we will 
discuss their expressive power, comparing them to the QTMs studied in 
the literature. The QTMs of B&V form a robust class, but meaningful 
computations are defined only for classical inputs (a single natural number 
with amplitude 1). Moreover, their QTMs “terminate” synchronously— 
either all paths in superpositions enter a final state at the same time, or all 
of them diverge. As a consequence, there is no chance to study—and give 
meaning—to infinite computations. More important, the class of “sensible” 
QTMs (in B&V’s terminology: the well-formed, well-behaved, normal form 
QTMs) is not recursively enumerable, since the constraint of “simultaneous 
termination” is undecidable. 

In Deutsch’s original proposal m, any quantum TM has an explicit ter¬ 
mination bit which the machine sets when entering a final configuration. 
While it is guaranteed that final probabilities are preserved, the observa¬ 
tion protocol requires that one checks termination at every step, since the 
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machine may well leave the final state (and change the tape). Deutsch’s ma¬ 
chines could in principle be used to define meaningful infinite computations, 
but we know of no such an attempt. 

In our analysis: (i) there is no termination bit: a quantum configuration 
is a genuine superposition of classical configurations; (ii) any computation 
path (that is, any element of the superposition) evolves independently from 
the others: any path may terminate at its own time, or may diverge; (iii) 
infinite computations are meaningful; (iv) we may observe the result of the 
computation in a way akin to Deutsch’s one, but with the guarantee that 
once a final state is entered, the machine will not change the correspond¬ 
ing “result” during a subsequent computation; (v) the class of QTMs is 
recursively enumerable, thus opening the door to a quantum computability 
theory which may follow some of the classical developments. 

2. Quantum Turing Machines 

In this section we define quantum Turing machines. We assume the reader 
be familiar with classical Turing machines (in case, see HD- 

In all the paper we shall assume that the tape alphabet E is finite and 
contains at least the symbols 1 and □: 1 will be used to code natural numbers 
in unary notation; □ will be the blank symbol. For any n G N, n will 
denote the string l n+1 . The greek letters cq/3, eventually indexed, will 
denote strings in E*; A will denote the empty string; af3 will denote the 
concatenation of a and (3. 

2.1. QTM. Given a tape alphabet E, we associate an extra tape symbol 
a to any a G E (including the blank □); E = {a | a G E} is the extra 
tape alphabet . Finally, E = E U E. As we discussed in the introduction, 
extra tape symbols will only appear in computations involving initial or final 
states. 

As in B&V, we assume that a machine move always implies a displacement 
of the machine head, to the left (L) or the right (R) on the tape; B = {L, R} 
is the set of the displacements. 

For / denumerable, £ 2 (/) is the Hilbert space of square summable, I- 
indexed sequences of complex numbers 

<£: J-HC| ^ \ cf )( C )\ 2 < oo 

Cei 

equipped with an inner product (. | .) and the euclidean norm ||<ft|| = y/ (4> | 0). 
See Appendix [A] for more details. 

Definition 1 (Quantum Turing machines). Given a finite set of states Q 
and an alphabet E, a Quantum Turing Machine (QTM) is a tuple 

M = (E, Q, Q s , Q t ,S,qi,q f ) 



where 
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• Q s C Q is the set of source states of M, and qi G Q s is a distinguished 
source state named the initial state of M; 

• Qt Q Q is the set of target states of M, and qj G Qt is a distinguished 
target state named the final state of M; 

• if we define Qo = Q \ ( Q s U Qt), then S = 5o U S s U St is the quantum 
transition function of M, defined as the union of the following three 
functions with disjoint domains: 

5o : ((Qo u Qs) xE)-> f((Qo U Qt) x E x D) 

S s '• (Qs x E) —> £ 2 (Q S x S x ID) 

S t : (Qt xE)4 f(Q t xExD) 


• the source transition function 5 S is defined by 


5 s (q s ,a)(q' s ,b,d) = 



if ( q' s ,b,d ) = ( q s ,a,R ) 
otherwise 


for every (q s , a) G Q s x X and every (<?', M)g2 s xSxD 
• the target transition function St is defined by 


St(qt,a)(q' t ,b, d) 


1 if (q' t ,b,d) = (q t ,a,R) 
0 otherwise 


for every (gt,a)GQxE and every G Qt xSxD 

• the main transition function 5o satisfies the following local unitary 
conditions 

(1) for any (q, a) <E (Qo U Q s ) x E 

X = 1 

(p,6,cf)E(QoUQt)xSxD 

(2) for any (g, a), (q r , a') <E (Q 0 UQ s )xE with (q, a) ^ (q 1 , a') 


X So(q',a')(p,b,d)*5o(q,a)(p,b,d) = 0 

(p,6,ci)e(QoUQt)xExD 

(3) for any (q, a, b), (q', a', b') e (Q 0 U Q s ) x E 2 


X $o(q', a')(p,b r , L)* 5 0 (q, a)(p,b, R) = 0 

pC(QoUQt) 

We remark that the domains and codomains of the three transition func¬ 
tions 5o, S s and St are disjoint. Indeed, if we define 

<$o — (Qo U Qs) x S S s = Q s x X <St = Qt x X 

7o — (Qo U Qt) x X T s — Q s x X 7t — Qt x X 

we can then write, in a compact way 

S x — $x £ 2 (T x x D) 


for x G {0, s, £} 
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and see that, for x, y G {0, s, t}, 

S x nS y = V) r x nT y = 0 

when x ^ y. 

2.2. Configurations. A configuration of a (classical) TM is a triple formed 
by the content of the tape, the state of the machine and the position of the 
tape head. As usual, we assume that only a finite portion of the tape contains 
non-blank symbols. We may therefore represent such a configuration as a 
triple (cqg,/3) G E* x Q x E* where: 

(1) q is the current state; 

(2) /? is the right content of the tape and its first symbol is the one under 
the head (we say also: in the current cell). That is, (3 = u/3', where 
the current symbol u is the content of the current cell (i.e. the one 
pointed by the tape head) and /?' is the longest string on the tape 
ending with a symbol different from □ and whose first symbol (if 
any) is written in the cell immediately to the right of the current 
cell; by convention, when the current cell and all the right content 
of the tape is empty, we shall also write (<a, g, A) instead of (<a, g, □); 

(3) a is the left content of the tape. That is, it is either the empty string 
A, or it is the longest string on the tape starting with a symbol differ¬ 
ent from □, and whose last symbol is written in the cell immediately 
to the left of the current cell. 

According to this definition, in a configuration (cqg,/3) the string a does 
not start with a blank □, and f3 does not end with a blank. In the following it 
will be useful to manipulate configurations which are extended with blanks 
to the right (of the right content) or to the left (of the left content). For 
this, we equate configurations up to the three equivalence relations induced 
by the following equations 

a ~i Da (3 (3n 

(<a, g, f3) (a ', g, (3') when a a! and (3 ~ r (3 r 

We now turn to QTMs. Observe first that while cells containing the blank 
symbol □ are considered empty, cells containing the extra symbol □ are not 
empty and should not be ignored on the left/right side of the tape. Moreover, 
in view of the particular evolution required for source/target states, and the 
special role of extra symbols, some of the triples (cqg,/3) cannot occur as 
actual configurations in a computation of a QTM. We limit our QTMs to 
configurations where extra symbols appear in af3 only when the current 
state is a source (target) state and, moreover all the extra symbols are 
immediately to the right (left) of the tape head. 

Definition 2 (configurations). Let M = (E, Q, Q s , Qt, S , go, Qf) be a QTM. 
A configuration of M is a triple (a, g, j3) G S* x Q x E* s.t. 

(1) if q 0 Q s U Qt, then af3 G £*, that is, the tape does not contain 
extra symbols; 
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(2) if q G Q s , then a G E* and /3 G E*E*; 

(3) dually, if q G Qt, then /3 G E* and a G E*E*. 

The set of configurations of M is denoted by Cm- A configuration of M is a 
source / target configuration when the corresponding state is a source/target 
state, moreover, it is a final/initial configuration when the current state is 
final/initial. By $m we shall denote the set of the final configurations of M. 

In the following, the index M in Cm and in the other names indexed by 
the machine might drop when clear from the context. 


2.3. Quantum configurations. As already discussed in the introduction, 
the evolution of a QTM is described by superpositions of configurations (as 
defined in Definition [2]). If B C Cm is a set of configurations, superpositions 
are elements of the Hilbert space £ 2 {B) (see, e.g., OUT]). Quantum config¬ 
urations of a QTM M are those elements of ^ 2 (Cm) with unitary norm. We 
remark that, since there is no bound on the size of the tape in a configura¬ 
tion, the Hilbert space of the configurations must be infinite dimensional. 


Definition 3 (quantum configurations). Let M be a QTM. The elements of 
the set qCM = {</> G ^ 2 (Cm) | J2ce£ M I^^OI 2 = 1} are the q-configurations 
(quantum configurations) of M. 

We shall use Dirac notation (see Appendix [A]) for the elements </>, ^ of 
qC m, writing them |</>), | ip). 


Definition 4. For any set of configurations B C Cm and any C G B let 
|C) : B -A C be the function 


|c> (D) = 


1 if C = D 

0 liC^D. 


The set CB(£>) of all such functions is a Hilbert basis for £ 2 {B) (see, e.g., jT5] ). 
In particular, following the literature on quantum computing, CB (Cm) is 
called the computational basis of £ 2 (C). Each element of the computational 
basis is called base q-configuration . 


With a little abuse of language we shall write |C) G \4>) when (/)(C) ^ 0. 
The span of CB(£>), denoted by span(CB(£>)), is the set of the finite linear 
combinations with complex coefficients of elements of CB(£>); span(£>) is 
not a Hilbert space, although £ 2 {B) is the (unique, up to isomorphism) 
completion of span(£>). Moreover, each unitary operator U on span(£>) has 
a unique unitary extension on £ 2 {B) 0. 

For a list of the main definitions, properties and notations on Hilbert 


presents a synoptic table of the so-called Dirac notation that we shall use in 
the paper. 


spaces with denumerable basis, see Appendix |A| In particular, subsection A.l 
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2.4. Time evolution operator. For any QTM M with alphabet X, and 
space of states Q, the step function 7 m • Cm xQxExD^ Cm is the 
map that, given a configuration of the tape and a triple (p, 6 , d) describing a 
“classical” step of a Turing machine, replaces the symbol in the current cell 
with the symbol p, moves the head on the d direction, and sets the machine 
into the new state p. Formally: 


7 M ((aw,q,u/3),p,v,d) ~ 


(awv,p, ft) 
(a,p, wv/3) 


when d — R 
when d — L. 


The evolution of a QTM M = (Q, X, S, go, Qf) can then be defined as a map 
on q-configurations. Following the three-parts definition of the transition 
function, let 


Cm = {(«, g, u/3) G Cm | (<?, r) G with x G {0, s, t}. 

It is easily seen that C^, C^, and C^ are a partition of Cm (they are 
pairwise disjoint, because <So, <$s, and St are pairwise disjoint). Therefore, 
given 

C = (<a, g, u/3) G Cm = 7 M (C,p,v,d) 

we can define 

Wm(|C)) = X u)(p,v,d) |C PiV>d ) when Ce£ x M . 

(p,v,d)eTx xB 


Proposition 5. Wm{ |C)) G span(CB(C)) ; for any C G C. Tden, VFm 
naturally extends to an automorphism on the linear space of q-configurations 

Wm • span(CB(C)) -A span(CB(C)). 


Proof. Let C and C p ^ v ^ be as in the definition of Wm- 

(1) Let C G <£? m . If (p, v) G 7o, then G C^ when p G Qt, and 

Cp^ v ^ G C^ otherwise. Thus, C P)V ^ G Cm, for every (p, x,d) G 

7o x D. 

(2) Let C G C|^, that is, u = a for some a G X and C = (a, g,a 7 / 3 ), 
where o,/3 ,7 G X*. If (p, x) G 7s, then d(g, d)(p, x, d) 7 ^ 0 iff x = a 
and d = i?, namely, Wm(|C)) = |C Pj a jR ) with <7^# = (aa,p, 7 / 3 ). 
Then, if the symbol a replaced by a was the last extra symbol on 
the tape, that is 7 = A, then (7 p?a ,-R G otherwise (7 p?a 

In any case, Wm(|C)) G span(CB(C)). 

(3) Let C G C^, that is, u — a for some a G X and C = ( 07 , g,/?), 

where cq/3 ,7 G X*. If (p, x) G 77, then d(g, a)(p, x, d) 7 ^ 0 iff x = a 

and d = R\ namely, W M (\C)) = |C p> a jR ) with = (a7a,p,/3). 

Then, G C^ and VFm(| c )) € span(CB(C)). 

Summing up, Wm( |C)) G span(CB(C)) in any case. Thus, Wm uniquely 
extends to an automorphism on span(CB(C)) by linearity: that is 

W M ( E MQ)= E k C W M ( |C)). 

Ce£ M Ce£ M 
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□ 

By completion, Wm extends in a unique way to an operator on the Hilbert 
space of q-configurations. 

Definition 6 (time evolution operator). The time evolution operator of M 
is the unique extension 

Um • ^ 2 (£m) ^ 2 (£m) 

of the linear operator Wm : span(CB(£)) -A span(CB(£)). 

Theorem 7. The time evolution operator of a QTM is unitary. 

Proof. The proof is a variant of the one given by B&V and by Nishimura and 
Ozawa [16j- In particular, we prove first that Um is an isometry of f? 2 (£), 
and then that, in this particular case, this implies that Um is unitary (which, 
in general, holds for finite dimensional Hilbert spaces only—in the infinite 
dimensional case an isometry might not be surjective). The full details of 
the proof are given in Appendix [Bj □ 

In Appendix |B] we shall not only show that the unitary local conditions 
imply the unitarity of the time evolution operator (i.e., Theorem [7|, but 
that they are also necessary. We remark that, this is not just a simple 
adaptation to our case of the already known proofs for B&V QTM; in fact, 
we also simplify the argument that allows to show that the isometry of Um 
implies its unitarity. 

Since the time evolution operator of a QTM is unitary, it preserves the 
norm of its argument, hence it maps q-configurations into q-configurations. 

Proposition 8. Let M be a QTM. If \<f>) G q(& then Um |0) £ q<& 

Definition 9 (initial and final configurations). A q-configuration \<f>) = 
e* |Q) is initial if all the Ci are initial, and it is final if all the Ci are final. 
Moreover, q^ is the set of final q-configurations, and we shall denote by |n) 
the initial configuration (A, qo,n). 

Definition 10 (computations). Let M be a QTM and let Um be its time 
evolution operator. For an initial q-configuration 10), the computation of 
M on |0) is the denumerable sequence {|0i)}zeN s.t. 

(1) \(po) = |0); 

(2) |0i) = t^|0). 

Clearly, any computation of a QTM M is univocally determined by its 
initial q-configuration. The computation of M on initial q-configuration 10) 
will be denoted by . 

Lemma 11. Let C G $ be a final configuration without extra symbols; in 
particular, let C = (/3,qf,a) with a(3 G and a = a\...ai. For every 
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k > 0, let 


C[k} = 


(fi,qf,ai . ..a k a k+ i ...a t ) 
(, 8,q f ,a 1 ...ain k ~ l ) 


for 0 < k < l 
otherwise 


that is, for k > 0, C[k\ is obtained from C by replacing the current symbols 
and the first k — 1 symbols to its right by the corresponding extra symbols. 

(1) |C[j]) = U J M l \C[\\), for every i,j > 0. More generally, for every 
i, j > 0 and every 0 G q£, 


(^TM.Ic[j]>) = W>|c[i]> 


(2) UjJ 1 |C[i]) G l 2 {ft \ $), for every i > 0 and j > 0. 

(3) Let 1 6) G l 2 (£ \ #). For every i , j > 0. we have that: 

(a) |C[i]> G W |0) only if i < j. 

(b) |C[ij) 0 U-i |0), that is, |0) G £ 2 (£\$) 


Proof 

(1) By the definition of the transition function, we can easily see that 

|C[i + 1]) = Um |C[i]), for every i > 0. Therefore, |C[j]) = U J M | C[0] ) = 
u m 1u m\ C 1°}) ~ Uif l \C\\]), for every i,j > 0. Then, since Um is 
unitoty, (ujj‘ 1 4.) , |Cp])) = W UlJ’ |C[j]) = | C[i]>. 

(2) Let D be any final configuration without extra symbols. By the 

previous item, |C[i]) , |D[k])^ = <C[i] | D[k + j + i]) = 0, for 

every k > 0, since j > 0 and then k + j + i > i. Therefore, 
since D[k\ 0 Uf^~ l |C[i] ), for any D and any k, we conclude that 
iy- |C[i])e< 2 (CU). 

(3) They are immediate consequences of the previous two items. Indeed, 
(u J M \(p) ,|C[i])) = (<f> | C[i - j]]) = Ofor i>j, and (ujj |</>), |C[i])) = 
(0 I C[i + j]]) = 0, for j > 0, since |0) G l 2 {£ \ 5). 

□ 


Remark 12. The previous lemma shows that the final configurations reached 
along a computation are stable and do not interfere with other branches 
of the computation in superposition, which may enter into a final con¬ 
figuration later. Indeed, given a configuration |0) = |0f) + |0 n f), where 
10f) G and |0 n f) does not contain any final configuration, let 0 = U l |0) = 
U l |0f) + U l |0nf). Any final configuration in U l |0 n f) contains less than i ex¬ 
tra symbols, while any final configuration C[k\ with k extra symbols in |0) is 
mapped into a configuration C[i+k\ with i+k extra symbols, without chang¬ 
ing the value associated to the configuration, since val(C[fc]) = val {C[k+i}). 
Moreover, C[k\ and C[i + k\ have the same coefficient in |0) and |0), respec¬ 
tively, since (0 | C[i + k]) = ( U l |0f), |C[i + k])) = (0 | C[k]). 
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2.5. A comparison with Bernstein and Vazirani’s QTMs: part 1. 

We refer to B&V for the precise definitions of the QTMs used in that paper. 
For the sake of readability, we informally recall the notion of what they call 
well formed, stationary, normal form QTMs (B&V-QTMs in the following). 

A B&V-QTM M = (E, Q,5,qo,qf) is defined as our QTM (with one 
source state and one target state only) with the following differences: 

(1) the set of configurations coincides with all possible classical config¬ 
urations, namely all the set E* x Q x £*. 

(2) no superposition is allowed in the initial q-configuration (it must be 
a classical configuration (a,q,/3) with amplitude 1); 

(3) let |C) be such an initial configuration and let 

k = min{j | U J M |C) contains a final configuration} 

If such a k exists, then (i) all the configurations in Um k |C) are final; 
(ii) for all i < k,UM % |C) does not contain any final configuration. 
We say in this case that the QTM halts in k steps in U^ |C); 

(4) if a QTM halts, then the tape head is on the start cell of the initial 
configuration; 

(5) there are no extra symbols and the transitions out of the final state 
or into the initial state are replaced by loops from the final state 
into the initial state, that is, S(qf,a)(qo,a,R) = 1 for every a G E. 
Therefore, because of the local unitary conditions, that must hold in 
the final state also, these are the only outgoing transitions from the 
final state and the only incoming state into the initial state, that is, 
5(qf,a)(q',a',d) = 0 if ( q',a',d ) ^ (q 0 ,a,R) and S(q', a')(qo, a, d) = 0 
if ( q',a',d ) ^ (q f ,a,R). 


Theorem 13. For any B&V-QTM M there is a QTM M f s.t. for each 
initial configuration |C) ; if M with input |C) halts in k steps in a final 
configuration |</>) = U\ |C), then U^,\C) = \(j>). 


Proof The QTM M f has the same states of M, only one source state, the 
initial state qo, and only one target state, the final state qf. Therefore, if 
M = (E, Q,S,q 0 ,q f ), we take M' (E, Q, {<?o}, {qf}^',q 0 ,q f ). 

The source part 5' and the target part S' t of the transition function S' of 
M' are uniquely determined by the definition of QTM. The function is 
instead the restriction of S to the domain <So — Q \ {Qf} x ^ that is, for 
each q ^ q/ and a G E, we have S' Q (q,a)(p,b,d) = S(q,a)(p,b,d ), for every 
(p,b,d) G 7o x D. The unitary local conditions hold for (%, since they hold 
for S and because, as already remarked, in a B&V-QTM, ii q ^ qf, then 
S(q, a)(qo, b, d) — 0, for every a, b G E and d G B. 

By construction, it is clear that for each i < k, s.t. \(j>\) = U l M , |C) is not 
final, then \(f>\) = U l M |C). □ 
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3. Quantum Computable Functions 

In this section we address the problem of defining the concept of quantum 
computable function in an “ideal” way, without taking into account any 
measurement protocol. The problem of the observation protocol will be 
addressed in Section [4j Here we show how each QTM naturally defines a 
computable function from the sphere of radius 1 in £ 2 to the set of (partial) 
probability distributions on the set of natural numbers. 

Definition 14 (Probability distributions). 

(1) A partial probability distribution (PPD) of natural numbers is a 

function P : N -A such that ( n ) — I* 

(2) If P( n ) = 1, P is a probability distribution (PD). 

(3) P and Pi denotes the sets of all the PPDs and PDs, respectively. 

(4) If the set {n : P(n) ^ 0} is finite, P is finite . 

(5) Let P', P" be two PPDs, we say that P' < P" (P' < P") iff for each 
n G N, P \n) < P "(n) (P \n) < P"(n)). 

(6) Let V — {PzjieN be a denumerable sequence of PPDs; V is monotone 
iff P^ < Pj, for each i < j. 

Remark 15. In the following, we shall also use the notation P(-L) = 1 — 
P( n )* By definition, 0 < P(-L) < 1, and a PPD is a PD iff P(_L) = 0. 

Since real numbers are a complete space, we have the following result. 

Proposition 16. Each PCP has a supremum, denoted by \_\V. 

Proof. For each n G N, the set V n — {P(n) : P G V} has a supremum 
|_| Vn- It is a trivial exercise to verify that (U V)(n) — \_\V n is indeed the 
supremum of V. □ 

We can now introduce the notion of limit of a sequence V = {P^}^. 

Definition 17. Let V — {Pi}i<N be a sequence of PPDs. If for each n G N 
there exists l n — hup-^ Pi(n), we say that lim^oo P^ = P, with P(n) = l n . 

The computed outputs of a QTM will be defined as the limit of the se¬ 
quence of partial probability distributions obtained along its computations. 

Definition 18 (probability and q-configurations). Given a configuration 
C = (a,#,/?), let val[C] be the number of 1 and 1 symbols in a(3. 

(1) To any q-configuration \(f>) — we associate the partial 

probability distribution s.t. P| ^(n) = Ece5,vai[C]=n l e c| 2 - 

(2) For any computation KH. = {l^i) • let P k m be the sequence of 

10) 

PPDs {P^j)}ieN- 

We now show the key property that a QTM computation yields monotone 
sequences of PPDs. In its simple proof we see at work all the constraints on 
the transition function of a QTM. First, once in a target state the machine 
can only change a (normal) symbol a into the corresponding extra symbol a; 
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as a consequence, the val of these configurations does not change. Second, 
when entering for the first time into a target state there are no extra symbols 
on the tape, for the form of the transition function 5 q, which is reflected 
in the constraints on the configurations (Definition [ 2 ]). Finally, in a final 
configuration the number of extra symbols counts the steps the machine 
performed since entering for the first time into a target state. These last two 
properties defuse quantum interference between final configurations reached 
in a different number of steps. 


Theorem 19 (monotonicity of computations). For any computation 
of a QTM M, the sequence of PPDs P k m is monotone. 


Proof. It is a direct consequence of the properties already remarked in 
Lemma 11 and Remark 12 Anyhow, let us see a direct proof in details. 

Let U be the time evolution operator of M. We prove that Vz Vn: 
P| 0 .)(n) < Pjj^.^n). Let us split any q-configuration in two parts, the 
sums of the final and of the non-final configurations: 


i0i)= e d c |c) + y> D |D) 


CeTi 




where Pi C $ are the final configurations in \(f)\) with non-null amplitude. 
By applying U we get 


U | 0 i) = E d cU |c> +J2 d ' B |B> + E d A |A> 

CeTi 


where J-[ C 5, and d' B / 0 for B E Jq, and 


U(Y J d D \D))= E 4 |B) + E d' A | A). 

Be T[ 

The sum on non-final configurations does not contribute to On the 

other hand, let ►F |C) = #{a : a G C}. For each C G Pp. 

(1) U |C) = |C') € S'; 

(2) val[C] = val[C"]; 

(3) *|C')=*|C) + 1>0; 

(4) for any other E € E% (i.e., E ^ C ), we have U |E) = E' ^ C'. 

As for the newly final configurations B 6 T[, we have |B) = 0. Therefore, 
none of the B E Ti is equal to any of the C' s.t. C') = U |C), and hence 


VU\4>;)= E \ d c\ 2 + E l^fll 2 

C&Ti B&T[ 

val [C]=n val[B]=n 


P IA>+ E l<4l 2 

BdT[ 

val[_B]=n 


□ 


We now turn to the definition of the computed output of a QTM compu¬ 
tation. The easy case is when a computation reaches a final q-configuration 
|f) G 5 (meaning that all the classical computations in superposition are 
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“terminated”)—in this case the computed output is the PD P|p). Of course 
the QTM will keep computing and transforming \fi) into other configura¬ 
tions, but these configurations all have the same PD. However, we want to 
give meaning also to “infinite” computations, which never reach a final q- 
configuration, yet producing some final configurations in the superpositions. 
In this case we define the computed output as the limit of the PPDs yielded 
by the computation. 

We need first the following lemma, whose proof is an easy consequence 
of the definition of limit for PPDs and of well known properties of limits of 
real-valued, non-decreasing sequences on natural numbers. 

Lemma 20. Let K = {|</>j)}^ be a monotone sequence of q-configurations, 
then the sequence of PPDs V = {P^}^ enjoys lim^ooPi^ = up- 

From the lemma and Theorem M we have: 

Corollary 21. Let K — {|0i)}ieN be a computation, then the sequence 
{P| 0 i)}ieN has limit lim^ooPi^ (denoted by lim K). 

The existence of the limit allows the following definition. 

Definition 22 (computed output of a QTM). The computed output of a 
QTM M on the initial q-configuration \<fi) is the PPD P = lim(nota¬ 
tion: M 10 ) -A P). 

Let us note that a QTM always has a computed output. 

Definition 23. Given a QTM M, a q-configuration | <jj) is finite if it is an 
element of span(CB(£|vi))- A computation = {|</>i)}zeN is finitary with 
computed output P if there exists a k s.t. is final and P|</> k ) = P. 

Proposition 24. Let = {|0i) be a finitary computation with com¬ 
puted output P. Then: 

(1) There exists a k, such that for each j > k, |0j) is final and P|^ = P; 

( 2 ) M \(p) p ; 

(3) P is a PD. 

Proof. Apply the definitions; P G Pi since |0|<) is final and has norm 1. □ 

With a computation several cases may thus happen: 

(1) K is finitary. The output of the computation is a PD, which is 
determined after a finite number of steps; 

(2) K is not finitary, and M^ -A P G Pi. The output is determined 
as a limit; 

(3) is not finitary, and -A P G P — Pi (the sum of the proba¬ 
bilities of observing natural numbers is p < 1). Not only the result is 
determined as a limit, but we cannot extract a PD from the output. 
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The first two cases above give rise to what Definition 25 calls a q-total 
function. Observe, however, that for an external observer, cases (2) and (3) 
are in general indistinguishable, since at any finite stage of the computation 
we may observe only a finite part of the computed output. 

For some examples of QTMs and their computed output, see Section [5] 


3.1. Quantum partial computable functions. We want our quantum 
computable functions to be defined over a natural extension of the natu¬ 
ral numbers. Recall that, for any n E N, n denotes the string l n+1 and 
that |n) = |(A,qo,n)). When using a QTM for computing a function, we 
stipulate that initial q-configurations are superpositions of initial classical 
configurations of the shape |n). Such q-configurations are naturally isomor¬ 
phic to the space £\ — {(f) : N -A C | J2 n eN l<K n )| 2 = 1} °f square summa- 
ble, denumerable sequences with unitary norm, under the bijective mapping 
iy(J2d k n k ) = X4 |l[k)- 

Definition 25 (partial quantum computable functions). 

(1) A function / : -A P is partial quantum computable (q-computable) 

if there exists a QTM M s.t. /(x) = P iff -A P. 

(2) A q-partial computable function / is quantum total (q-total) if for 
each x, /(x) E Pi. 

QCT is the class of partial quantum computable functions . 


4. Observables 

While the evolution of a closed quantum system (e.g., a QTM) is re¬ 
versible and deterministic once its evolution operator is known, a (global) 
measurement of a q-configuration is an irreversible process, which causes 
the collapses of the quantum state to a new state with a certain probability. 
Technically, a measurement corresponds to a projection on a subspace of the 
Hilbert space of quantum states. For the sake of simplicity, in the case of 
QTMs, let us restrict to measurements observing if a configuration belongs 
to the subspace described by some set of configurations £>. The effect of 
such a measurement is summarised by the following: 

Measurement postulate 

Given a set of configurations B C (£, a measurement observ¬ 
ing if a quantum configuration \cf)) = J2ce£ ec IQ belongs 
to the subspace generated by CB (£>) gives a positive answer 
with a probability p — J^ceB \ e c\ 2 i equal to the square of 
the norm of the projection of |</>) onto £ 2 (£>), causing at the 
same time a collapse of the configuration into the normalised 
projection ^2ceBP~ le C |Q; dually, it gives a negative answer 
with probability 1 — p — J2c^B \ e c\ 2 and a collapse onto the 
subspace £ 2 ((£ \ B) orthonormal to £ 2 (B), that is, into the 
normalised configuration J2c^b(^ ~ P)~ le C IQ- 
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Because of the irreversible modification produced by any measurement 
on the current configuration, and therefore on the rest of the computation, 
we must deal with the problem of how to read the result of a computation. 
In other words, we need to establish some protocol to observe when a QTM 
has eventually reached a final configuration, and to read the corresponding 
result. 

4.1. The approach of Bernstein and Vazirani. We already discussed 
how B&V’s “sensible” QTMs are machines where all the computations in 
superposition are in some sense terminating, and reach the final state at the 
same time (are “stationary”, in their terminology). More precisely, Defini¬ 
tion 3.11 of B&V reads: ”A final configuration of a QTM is any configura¬ 
tion in [final] state. If when QTM M is run with input x, at time T the 
superposition contains only final configurations, and at any time less than T 
the superposition contains no final configuration, then M halts with running 
time T on input x. ” 

This is a good definition for a theory of computational complexity (where 
the problems are classic, and the inputs of QTMs are always classic) but it is 
of little use for developing a theory of effective quantum functions. Indeed, 
inputs of a B&V-QTM must be classical—we cannot extend by linearity a 
B&V-QTM on inputs in since there is no guarantee whatsoever that on 
different inputs the same QTM halts with the same running time. 

4.2. The approach of Deutsch. Deutsch m assumes that QTMs are 
enriched with a termination bit T. At the beginning of a computation, T is 
set to 0, and the machine sets this termination bit to 1 when it enters into a 
final configuration. If we write |T = i) for the function that returns 1 when 
the termination bit is set to i, and 0 otherwise, a generic q-configuration of 
a Deutsch’s QTM can be written as 


\d>) = |T = 0) <g> ]T e c |C) + |T = 1) <g> Y, d D |D> 

Des 

The observer periodically measures T in a non destructive way (that is, 
without modifying the rest of the state of the machine). 

(1) If the result of the measurement of T gives the value 0, | <f>) collapses 

(with a probability equal to l e c| 2 ) to the q-configuration 

T = 0) ® Xc£g e c IQ 

l e d 2 

and the computation continues with 

(2) If the result of the measurement of T gives the value 1, |</>) collapses 

(with probability \^d\ 2 ) to 

| ,//\ |T = 1) ® l D ) 
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and, immediately after the collapse, the observer makes a further 
measurement of the component 
a final configuration. 


,d D D 


in order to read-back 


Note that Deutsch’s protocol (in a irreversible way) spoils at each step the 
superposition of configurations. The main point of Deutsch’s approach is 
that a measurement must be performed immediately after some computation 
enters into a final state. In fact, since at the following step the evolution 
might lead the machine to exit the final state modifying the content of 
the tape, we would not be able to measure at all this output. In other 
words, either the termination bit acts as a trigger that forces a measurement 
each time it is set, or we perform a measurement after each step of the 
computation. 


4.3. Our approach. The measurement of the output computed by our 
QTMs can be performed following a variant of Deutsch’s approach. Because 
of the particular structure of the transition function of our QTMs, we shall 
see that we do not need any additional termination bit, that a measurement 
can be performed at any moment of the computation, and that indeed we 
can perform several measurements at distinct points of the computation 
without altering the result (in terms of the probabilistic distribution of the 
observed output). 

Given a q-configuration |Q = |0f) + |0 n f), where |0f) G £ 2 ($) and \4> n f) E 
t? 2 (£ \ 3% our output measurement tries to get an output value from 1 (j>) by 
the following procedure: 

(1) first of all, we observe the final states of |</>), forcing the q-configuration 
to collapse either into the final q-configuration |0f) / |||0f)||, or into 
the q-configuration \(f> n f) / |||</>nf)|| 5 which does not contain any final 
configuration; 

(2) then, if the q-configuration collapses into \<fif) /|||</>f)||, we observe 
one of these configurations, say |C), which gives us the observed 
output val[C] = n, forcing the q-configuration to collapse into the 
final base q-configuration (e c /|e c |) |C); 

(3) otherwise, we leave unchanged the q-configuration \(f> n f) / |||0 n f)|| ob¬ 
tained after the first observation, and we say that we have observed 
the special value _L. 

Summing up, an output measurement of |</>) may lead to observe an output 
value n G N associated to a collapse into a base final configuration |C) G | <j>) 
s.t. val[</>] = n or to observe the special value _L associated to a collapse into 
a q-configuration which does not contain any final configuration. 

Definition 26 (output observation). An output observation with collapsed 
q-configuration |Q and observed output x G N U {T} is the result of an 
output measurement of the q-configuration |</>) = J2ce<t ec IQ- Therefore, 
it is a triple |</>) |Q s.t. 
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(1) either x — n G N, and 


= “j ~r |C) with |C) G |0f) and val[C] = n 


(2) or x = i, and 
i \M 


where \4> nf ) = ^ e c |C) and |||0 n f)|| + 0 

cn 


Ill0nf)|| 

The probability of an output observation is defined by 

( \ec\ 2 if x = n G N 


Pr{| 4>) lx IV’)} = 


l 


II l*/*nf) lr 


if x = _L 


Remark 27. Let e|C) ^ 10), with C G 5 and val[C] = n. By definition, 
x — n and \(j>) = (e/|e|) |C); moreover, Pr{e |C) ^ |0)} = |e| 2 . 


Remark 28. For every distinct pair of output observations |0) ^ |i/q) and 
10) 0^2 1020 we have that and ^2 are in the orthonormal subspaces 
generated by the two disjoint sets £> i ,£>2 50 £, where — {C G £ | |C) G 

IV’i)}- 


Definition 29 (observed run). Let M be a QTM and Um its time evolution 
operator. For any monotone increasing function r : N -A N (that is, r(i) < 
r(j ) for i < j): 

(1) a r-observed run of M on the initial q-configuration |</>) is a sequence 
{|0i)}ieN s.t.: 

(a) |0 O ) = 10); 

(b) Um 10h) ixi |0h+i), when h = r{i) for some i G N; 

(c) |0h+i) = U M 10h) otherwise. 

(2) A finite r-observed run of length k is any finite prefix of length k 
of some r- observed run. Notation: if R — {|0i)}ieN, then R[k] = 
{|0i) }i<k- 


Remark 30. We stress that, given an R = {|0i)}ieN : 

(1) either it never obtains a value n G N as the result of an output 
observation, and then it never reaches a final configuration; 

(2) or it eventually obtains such a value collapsing the q-configuration 
into a base final configuration u |C) s.t. \u\ = 1 and val[C] = n, and 
from that point onward all the configurations of the run are base 
final configurations u |Cj) = uU J |C) s.t. val[Cj] = n, and all the 
following observed outputs are equal to n (see Remark 27). 


Definition 31. Let R — { |0i)be a r-observed run. 

(1) The sequence s.t. |0h) i Xi |0h+i)? with h — r(i), is the output 

sequence of the r-observed run R. 

(2) The observed output of R is the value xgNU { _L} (notation: R 
defined by: 
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(a) x = n G N, if X{ — n for some i G N; 

(b) x = JL otherwise. 

(3) For any fc, the output sequence of the finite r -observed run R[r{k)], 
is the finite sequence { Xi } i<T (&) and x & is its observed output. 


Definition 32 (probability of a run). Let R = {|0i)}ieN be a r -observed 
run. 


(i) 


For fceN, the probability of the finite r-observed run R[k ] is induc¬ 
tively defined by 
(a) Pr{i?[0]} = 1; 


(b) Pr{i?[fc+1]} 


Pr{i?[fc]}Pr{|0 k )U |4+i)} 


when k = r(i) 
for some i G N 


iPr {R[k}} 

(2) Pr {R} = Unifc^oo Pr {R[k]}. 


otherwise 


We remark that Pr{i?} is well-defined, since 1 > Pr{R[i]} > Pr{R[j]} > 0, 
for every i < j. Therefore, 

1 > Pr{i?} = lim Pr{R[k}} = inf{Pr{i?[/c]}}^ G ^ > 0. 

k—> oo 


Remark 33. Let R = {0z}zeN be a r -observed run s.t. R for some n G N. 
As observed in Remark |30| for some k, we have i?[r(fc)] and R[r(k) + 1] 
moreover, for i > k, |0j) = n|Cj) with |n| = 1, Ci G 5, and val[Cz] = n. 
As a consequence, P r{R[k + 1]} = P r{R[i]} = Pr{i?}, for i > k (since, by 

10r(i)} -i-n |0r(i)+l)} !)• 


Remark 


27 


Definition 34 (observed computation). The r-observed computation of a 
QTM M on the initial q-configuration 10), is the set /C^ r of the r-observed 
runs of M on |0) with the measure Pr : V(IC r ) -A C defined by 


Pr£= Pr{R} 
ReB 

for every B C JCM. . 


By /C[fc]j^ r we shall denote the set of the finite r-observed runs of length 
k of M on |0), with the measure Pr on its subsets (see Definition 34). 

It is immediate to observe that the set /Cj^ r naturally defines an infinite 
tree labelled with q-configurations where each infinite path starting from 


the root |0) correspond to a r-observed run in K M 


I0)r* 


Lemma 35. Given R\,R 2 G /Cj^ jT , with Ri = {0i,z}zeN {</>2,z}zeN = R 2 , 
there is k > 0 s.t 

(1) 01 = 02,i for i < r{k), that is, Ri[r(k)] = 7?2[t(/c)]; 

( 2 ) for i > r{k), the q-configurations 0 7 ^ 02 ,z are in too orthonormal 
subspaces generated by two distinct subsets of £. 
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Proof. Let R\[h\ — R 2 [h\ be the longest common prefix of R\ and i? 2 - Since 
they both starts with |</>), such prefix is not empty; moreover, by the defini¬ 
tion of r-observed run, it is readily seen that h = r(fc), for some k. By con¬ 
struction, </>i fh +i / 02 ,/i+i and (f) h f Xj 1 , for j = 1,2, with <f h = (j) hh = 
,h- Moreover, at least one of the two q-configurations l^ph+i) , 1^2,h+i) is 
a final base q-configuration u |Ci); for instance, let |^i 5 h+i) = Ri |Ci). 

Let us take B a ,i = {C G £ | |C) G |} 5 for a — 1,2 and i > h. We 
prove that B\^ H £> 2 ,i = 0, for i > h. First of all, this holds for z = h + 1, by 
Remark [28| Then, we distinguish two cases: 

(1) 1^2,h+i) = ^2 |C 2 ) is a final base configuration. We have that l^i^+i) = 
u\ U l ~ h |Ci) 7 ^ 1/2 U l —h |C 2 ) = |^ 2 ,i+i), for a — 1,2 and i > h (by 
Remark [30] and the fact that Um is injective, since we are already 
remarked that C\ 7 ^ C 2 ). 

(2) 1^2,h+i) £ ^ 2 (£\30. Every |C) G |^ 2 ,h+i) contains less than i — h 
extra symbols, while I'i/’pi+i) = u\ U l ~ h |Ci) contains at least i — h 
extra symbols. Therefore, B\^+\ H £> 2 , 2+1 = {U l ~ h |Ci)} H £> 2 , 2+1 = 0. 

□ 


Lemma 36. = {<^}^ be the computation of the QTM M on 

the initial q-configuration \(f>) and 1C r r-observed computation on the 
same initial configuration. For every k G N, we have that 


1+) = E Pr W I^R> 


ReK[k] 


M 

1 4>),r 


where |^r) is the last q-configuration of the finite run R of length k. 


Proof. By definition, <f = </>o and R = {</>} with Pr{i?} = 1 and ifn — </>, is 
the only run of length 0 in /Cj^ r . Therefore, the assertion trivially holds for 
k = 0 . 

Let us then prove the assertion by induction on k. By definition and the 
induction hypothesis 

\A+i) = U M 14} = E U M |4> 

We have two possibilities: 

(1) k 7 ^ r{i) for any i. In this case, there is a bijection between the runs 
of length k and those of length k + 1, since each run R' G IC[k+ 1]|^ r 

is obtained from a path R G JC[k\ T with last q-configuration |^r), 
by appending to R the q-configuration |^r/) = Um I^r)- Moreover, 
since by definition, P r{R'} = Pr{R}, we can conclude that 

I0k+1> = E Pr W Um IV’r) = E PrWl<M 

Rem^, T R'eK[k+ i]+ T 
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(2) k = r(z), for some i. In this case, every R G JC[k\^ r with last q- 
configuration |^r) generates a run R' of length k + 1 for every output 
observation |^r) |^r/), where R' is obtained by appending |^r/) 

to R. Therefore, let R = {\^\)}i<k £ /C[fc]j^ T and 

= 1 | iV’k) iV’k+l)} 

by applying Definition [26j we easily check that 

Um I^r) = E Pl'd'i/’R) Irr IV’R')} IV’R') 

R'eB R 

Thus, by substitution, and Pr{i?} Pr{|^R) 4-cc |^r')} = Pi'{-R / } 

lA+l) = E ^{R}u M m 

R em^, T 

= E E Pr{,R} Pr{|^R) IUr')} IUrO 

*«:[*]" , r «'6B» 

E Pr{^}|^R-> 

U R£ic[k\M T 

= e i*'> 

R'e/c[/c+i]^ ;r 

since UfleJC[fc]" >iT &R = £[* + 1]JJ )T . 

□ 


We are finally in the position to prove that our observation protocol is 
compatible with the probability distributions that we defined as computed 
output of a QTM computation. 


Theorem 37. Let Kfff = {0z]teN be the computation of the QTM M on 


M) 


the initial q-configuration |0) and 1C ^ r the r-observed computation on the 
same initial configuration. For every n G N: 

(1) P|0 k )(n) = Pr{i? € /C[fc]|^ iT | R In}, for every k = r(i), with i € N; 

( 2 ) P K M(n) = Pr{RelC$ )>T )Rl n }. 


Proof. By Lemma 36, we know that | <f>\f) = P r {^} I^r), where 

AJ|0) jT 


|^r) is the last q-configuration of R. Since k — r(z), for some i, we also 
know that either fjR G £ \ ^ or |^r) = ur |Cr) with Cr G $ and \ur\ = 1 . 
Therefore, 


E y/Pr{R}u R \C R ) 

ReB[k,n] 


p |4)U) 


2 
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where 

B[k,n] = {RelC[k}^ T \RU} 

= {Re K\k]\4,),T I IV’r) = u R |C R ) with val [Cj*] = n} 

By Lemma [35j we know that for every i?i, i?2 £ £>[&, n], we have l C Ri) ^ 
|Cr 2 ). Therefore 

P|4)H= E P r{R}\u R \ 2 = E Pr {G= Pr<B[fc,n] 

ReB[k,n\ ReB[k,n] 

since \ur\ = 1 . Which concludes the proof of the first item of the assertion. 

In order to prove the second item, let B[u,n\ = {R G /C^ r | i? ^ n }. We 
have that 


Pr{i? € K M 


l<W, T 


R In} = E P'TG = , lim E 

zJ k^ 00 z ' 

Re£>[a;,Ti] 


Pr{i?} 


-ReS[w,n] 

R[r(k)]U 


^ = lim 

/ c — t >-00 




RE£>[ct;,n] 

H[r(fe)]j n 


( **) = lim 

k—>oo 


E 


Prli? 7 } = lim Pr B[r(k), n] 

k—> 00 


il'GiS[r(A;),n] 

since: (*) Pr{i?} = Pr{i?[r(fc)]}, when R[r(k)] i n (see Remark [33]) ; (**) 
there is a bijection between the sets S[r(k),n] and {R G B[uj, n] \ R[r(k)] i n } 
(see Remark |30|) mapping every R' G B[r{k)^n\ with last q-configuration 
u |C) into R — {ipR : i}ie N G {i? G £>[o;,n] | i?[r(fc)] s.t R' = R[r(k)] and 

\^ R;i ) = u U^ T{k) \C) for i>r(k). 

Therefore, by the (already proved) first item of the assertion 


Pr { R G K \}),r \Rln}= hm Pr B[r(k),n} = Urn P|* T(k) )(«) = P K * H 


□ 


5 . Remarks on the expressive power 

5.1. Computable configurations. Since QTMs represent (ideal) physi¬ 
cally realisable devices, we should constrain the complex numbers in the 
time evolution operator to be computable (see, e.g., Remark 9.2 in mi for 
a discussion). 

Definition 38 (computable numbers). A real number x is computable if 
there exists a deterministic Turing machine that on input n computes a 
binary representation of an integer m G Z such that |pr — x\ < 

The computable complex numbers C are those complexes whose real and 
imaginary parts are both computable. 


TOWARDS A THEORY OF QUANTUM COMPUTABILITY 


23 


Definition 39 (computable QTM). A QTM is computable iff for any (g, a) G 
<So and every (p, 5, d) G 7o x B, we have 5o(g, a)(p, 5, d) G C. 

Observe that, being (Qo U Qt) xExD finite, in a computable QTM the 
element 5o (g, a) G i 2 ((Qo U Qt) xExD) may be “effectively presented” in 
an obvious way. 

We postpone to a subsequent paper a full treatment of computable QTMs, 
and especially of the computability theory they may engender. We make 
here only some simple, preliminary remarks. 

First, computable QTMs form a recursive enumerable class. Indeed, any 
complex number e G C may be described by the index of the (classical) TM 
computing it (write ffefl G N for this). Moreover, for any (q,a) G <So we 
have a classical TM enumerating So(q,a) (that is, producing the family of 
the indexes f[e~[] of the TMs computing the amplitudes). 

Second, the time evolution operator Um of a computable QTM defines 
a classically computable function on the (code of) quantum configurations. 
We spell this out in case of finite q-configuration0 

Let us extend the alphabet S of the QTM M into E c = E U Q U {(,) *}. 
Any finite q-configuration | <j>) = Y^i=i e i \ ( a 'n Q'j P\)) ma y be coded by the 
string 

[[10)11 = ^ [[eilK aqg 1 fifi) * • • • ★ [ ]j^(a n q n f3 n )* G £*. 

Let us denote by [[qC/ull the set of such codes of finite q-configurations. 
The function : ffqCjwfH ([q^Affl defined by: 

n 

^ |(«i,q',/?i)))H 

2=1 

is intuitively computable, and therefore, by Church’s thesis, is computable. 
Therefore it is only a matter of routine to prove the following theorem: 


Theorem 40 (classical soundness). Let M be a computable QTM such that 
for each finite input \<f) the corresponding computation is finitary. There is 
a classical partial computable function 


s.t. for each finite 

Comp M ([[|</>)]l) = 


Comp M : [fqCji/P [fq^AfU 
G qC, M -A P iff there is 


s.t. P = P 


W 


and 


5.2. A comparison with Bernstein and Vazirani’s QTMs: part 2. 


In view of Theorem 13, we may say that our QTMs generalise B&V-QTMs, 
which may be simulated. The general framework, however, is substantially 
modified and the “same” machine behaves in different ways in the two ap¬ 
proaches. We give two simple examples of this, before concluding the paper. 


4 The argument may be generalised to infinite q-configurations—that is, q- 
configurations in which there is an infinite number of non-zero configurations in 
superpositions—provided these infinite superpositions are recursively enumerable. 
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Figure 1. The transition function S 

In this section we shall use a pictorial representation of QTMs, via a graph 
for the transition function S. When 5(q, a)(p, b,d) = x with x / 0 we draw 
the labelled arc as in Figure [I] 





Figure 2. Reversible TM a la Bernstein and Vazirani 


Example 41 (classical reversible TM with quantum behaviour). Let 

M be the reversible TM represented in Figure [2j where a E {□, 1}. M is a 
B&V-QTM indeed. 

If we feed M with a non classical input, e.g. |^) = ^ ID + ^ |3), then 
M fails to give an answer according to B&V’s framework, since B&V-QTMs 
always presuppose a classical input. 

If we transform M in our formalism (see Theorem 13), we obtain the 
QTM in Figure [3| From the definition of computed output, we have that 
M 


v^li)+v^l4) 


{\ : 2}; namely, with probability \ the QTM halts with 


output 2; while with probability \ it diverges. 


Example 42 (A PD obtained as a limit). The following example shows a 
machine which produces a PD only as an infinite limit. Let us consider the 
QTM in Figure[4j where £ = {$, 1, □}, a E {1, □}, p is a target state, and s 
is a source state. A simple calculation show that -A {1 : n+1}; namely, 
the machine M on input n produces with probability 1 the successor n + 1. 
We can see also that the PD {1 : n+1} is obtained only as a limit. Of 
course we must not wait an infinite time to readback the result! A correct 
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Figure 3. Reversible TM 

way to interpret this fact, is that for each n G N, each e G (0, there exist 
a natural number j s.t. P[/j| n ) (n + 1) > 1 — e. 


a~a,R,1A/2 $~$,R,1 



6. Conclusions and further work 

We find surprising that in the thirty years since m a theory of quantum 
computable functions did not develop, and that the main interest remained 
in QTMs as computing devices for classical problems/functions. This in 
sharp contrast with the original (Feynman’s and Deutsch’s) aim to have a 
better computing simulation of the physical world. 
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As always in these foundational studies, we had to go back to the ba¬ 
sics, and look for a notion of QTM general enough to encompass previous 
approaches (for instance, simulation of B&V-QTMs, Theorem 13), and still 
sufficiently constrained to allow for a neat mathematical framework (for in¬ 
stance, monotonicity of quantum computations, Theorem 19, a consequence 
of the particular way final states are treated in order to defuse quantum 
interference once such states are entered). While several details of the pro¬ 
posed approach may well change during further study, we are particularly 
happy to have a recursive enumerable class of QTMs. This may allow a 
fresh look to the problem of a quantum universal machine, and, therefore, 
to obtain some of the “standard” theorems of classical computability theory 
(s-m-n, normal form, recursion, etc.). These themes, as well those related 
to the various degrees of partiality of quantum computable functions (see 
the brief discussion after Proposition [24]) will be the subject of forthcoming 
papers. 
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Appendix A. Hilbert spaces with denumerable basis 


Definition 43 (Hilbert space of configurations). Given a denumerable set 
B , with £ 2 {B) we shall denote the infinite dimensional Hilbert space defined 
as follow. 

The set of vectors in £ 2 {B) is the set 

L|^:£^C,^|^(C)| 2 <ool 

l ces J 

and equipped with: 

(1) An inner sum + : l 2 (B) x £?{B) -A £ 2 (B) 
defined by (0 + 0 )(C ) = 0(C) + 0(C); 

(2) A multiplication by a scalar • : C x £ 2 (B) -A l 2 (B) 
defined by (a • 0)(C) = a • (0(C)); 

(3) An inner product|j (•, •) : 1 2 {B) x £ 2 {B) -A C 
defined by (0,0) = J^ceB 0( c ')*0( c '); 

(4) The Euclidian norm is defined as ||0|| = (0,0). 

The Hilbert space £ 2 = f(N) is the standard Hilbert space of denu¬ 
merable dimension—all the Hilbert spaces with denumerable dimension are 
isomorphic to it. £ 2 is the set of the vectors of £ 2 with unitary norm. 


Definition 44 (computational basis). The set of functions 


CB(£>) = {|C) :CeB, |C) : B ^ C} 


such that for each C 


|c> (D) = 


1 if C = D 
0 if C ± D 


is called computational basis of £ 2 {B). 


5 The condition ^2 CeB 0(C) | 2 < oo implies that ^ CeB 0(C)*0(C) converges for every 
pair of vectors. 
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We can prove that m- 

Theorem 45. The set CB (B) is an Hilbert basis of £ 2 (B). 

Let us note that the inner product space span(CB(£>)) defined by: 

span(CB(£>)) = | a G C,# G CB(£>), n G N 

L=1 

is a proper inner product subspace of £ 2 (B), but it is not an Hilbert Space 
(this means that CB (13) is not an Hamel basis of £ 2 (B)). 

The completion of span(CB(£>)) is a space isomorphic to £ 2 (B). 

By means of a standard result in functional analysis we have: 

Theorem 46. 

(1) span(CB(£>)) is a dense subspace of £ 2 (B); 

(2) £ 2 (B) is the (unique! up to isomorphism) completion o/span(CB(£>)). 

Definition 47. Let V be a complex inner product space, a linear application 
U : V -A V is called an isometry if ( Ux,Uy) = (x,y), for each x,y G V; 
moreover if U is also surjective, then it is called unitary . 

Since an isometry is injective, a unitary operator is invertible, and more¬ 
over, its inverse is also unitary. 

Definition 48. Let V be a complex inner product vectorial space, a linear 
application L : V -A V is called bounded if 3c > 0 Mx \Lx\ < c\\x\\. 

Theorem 49. Let V be a complex inner product vectorial space, for each 
bounded application U : V -A V there is one and only one bounded application 
U*:V^Vs.t. (x,Uy) = (U*x,y ). We say that U* is the adjoint ofU. 

It is easy to show that if U is a bounded application, then U is unitary 
iff U is invertible and (7* = U~ l . 

Theorem 50. Each unitary operator U in span(CB(£>)) has an unique ex¬ 
tension in £ 2 (B) 0. 

A.l. Dirac notation. We conclude this brief digest on Hilbert spaces, by 
a synopsis of the so-called Dirac notation, extensively used in the paper. 


mathematical notion 

Dirac notation 

inner product (0, ip) 

(<p i $) 

vector (p 

i*> 

dual of vector cp 
i.e., the linear application d^ 
defined as d^'ijf) — ((f),'ip) 

(01 

note that (cp \ ip) = (<p>\ ( ^)) 


Let L be a linear application, with (cf)\ L \^>) we denote {(/) \ Lif). 
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Appendix B. Proof of Theorem 0 

In this section we shall give the proof of Theorem [7] in full details. Follow¬ 
ing Bernstein and Vazirani |3j and Nishimura and Ozawa m, we shall prove 
a stronger result indeed, that the local unitary conditions are not only suf¬ 
ficient to obtain a unitary time evolution, but also necessary (Theorem |60|) . 

B.l. Pre-QTM. A pre-QTM is a tuple M = (£, Q, Q s , Q t , S, qi, qf) for 
which all the requirements demanded for a QTM (Definition |T]) but the 
local unitary conditions hold. In other words, while 5 S and St are defined 
and constrained as for QTM’s, we do not have any condition on Sq. Ground 
and quantum configurations, step function and time evolution operator of 
pre-QTM’s are defined as for QTM’s. In the following, M will denote a 
pre-QTM. 

B.2. Basic results. Let us start with some basic results which hold in any 
Hilbert space. As usual, if U is an operator, U* denotes its adjoint. 

Lemma 51. Let U : ^ 2 (£m) —^ ^ 2 (Cm)- 

(1) U is an isometry iff it is left-invertible and its adjoint is its left- 
inverse, that is, U*U = 1. 

(2) If U is an isometry, then P — UU * is an orthonormal projection 
and 

(a) ||P|0)|| < |||0)||, for every cf; 

(b) ||P |0)|| = |||0)|| tfP|0) = |0); 

(c) P=1 iff ||P |C) || = <C| P |C) = 1, for every C& <t M ■ 

Proof. 

(1) (£/J0 ), U 10)) = (U*U 10), 10)) = (010) iff U*U = 1* = 1. 

(2) P is a projection when PP — P and, moreover, it is orthonormal 

when it is self-adjoint, that is P* = P. Equivalently, P is a projection 
iff, for every 0, we have \(fj) = + \<fio) for some (unique) <f>\ and </>o 

s.t. P \(f>i) = (f>\ and P \(f>o) = 0; it is orthonormal when (</>i | </>o) = 0 
indeed. 

P — UU * is clearly self-adjoint. Moreover, PP = UU*UU* = 
UU * = P, since U is an isometry and then U*U = 1. Thus, P is an 
orthonormal projection. 

Let | (jj) — |0i) + \<fio) be an orthonormal decomposition as above. 

|||0)|| = 1110011+23? (0! | 0o) + |||0o)|| = ||P|0l)|| + |||0o)|| = IIP 10)11 + 
IH0o>||. 

(a) ||P|0)|| < ||P|0)|| + |||0o)|| = 1110)11 

(b) \\P 10)j| = jj 10)|| iff || |0o)II = 0, that is, iff |0 O ) = 0. 

(c) By the previous item, P |0) = |0) for every 0, iff ||P |0) || = |||0)|| 
for every 0, that is, iff ||P |C)|| = |||C)|| = 1 for every C G 
that is, iff ||P|C)|| = (P|C),P|C)) = (C|P|C) = 1, for every 
C G Cm- 

□ 
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B.3. Reverse transitions. The reverse step function of M is defined by 


7e ,Q({av R ,p,wv L /3),q,u,d) ~ 


(av R ,q,uf3) 
(a, q, uwvlP) 


when d — L 
when d — R 


Let us say that an R or L step of 7 is an iLreverse or L-reverse step. We 
see that an R/L -reverse step of 7 revert an i?/L step of the step function 
7 . While in 7 both the L-step and the R -step replace the same symbol, 
the current symbol of the configuration, in 7 the symbols replaced by the 
Lbreverse step and by the L-reverse step are in different positions. We have 
then a current L-reverse symbol vr and a current R -reverse symbol vr. 


Lemma 52. Let C ~ (a\Wi,q,uPi) and D ~ (a 2 VR,p, w 2 vrP 2 ). 

7 (D,q,u,d) = C iff 7 (C,p,v d ,d) = D 

Proof. When d — R, we have C ~ ( 02 , g, uw 2 vrP 2 ) — l(P,q,u,R) iff 
aiwi ~i a 2 and w 2 v L p 2 -r Pi iff 7 (C,P,vr,R) ^ {ol\WiVr, q, pi) ~ D. 
When d — L, we have C — (c£ 2 vrw 2 , q,up 2 ) ~ j(D,q,u,L) iff a\ ^7 a 2 vr 
and w\ — w 2 and p 2 ~ r f3\ iff 7 (C,p, vr, L) ( 07 , q , w\VrPi) — D. □ 

Lemma 53. Let C[q,u\ be the configuration obtained by substituting the 
state q and the symbol u for the current state and the current symbol of the 
configuration C, and let C^dl ~ 7 ( 7 ^, p, v, d\), q,u, d 2 ). 

(1) If d = di = d 2 , then C[q,u} = C p q ^ d d , for every (p,v) eQxS. 

(2) When d\ d 2 , there is Cd 1 d 2 [l v \> which does not depend on p, 
s-t. C dl d 2 [q,u,v] = <77 *for every p <E Q. 

Moreover, let (q f , v! , v'\ (q" , u", v") GQxSxS and d\,d 2 , d! x , d 2 E B^ with 
d\ 7 ^ g ?2 and d[ d 2 . If the tape of the configuration C contains at least a 
non-empty cell in addition to the current one (that is, C (A ,q,u\)), then 

(3) CV, u'] = C[q", u"] iff (q ', u') = (g", u")/ 

(4) C[g',u'j / C'd 1 d 2 b /, ,w ,/ ,v"]; 

(5) <V 1 d 2 [g / ,’U / ,'u / ] = C' d / i ^[g",«",v"] iff(di,d 2 ) = K,4) and (q 1 , u', v') = 
W■> u " 1 v") 

Proof. Let C ~ (aziwi,q,uw r z r f5). By computing <7’“’^, we see that (1) 
and ( 2 ) hold with 

C[q',u'] ~ (aziwi,q',u'w r z r /3) 

ClrW,u',v\ ~ (a, </, u'wivw r z r (3) 

CRL[q',u',v\ ~ (aziwivw r ,q',u'(3) 

From which, we can prove the following items. 

(3) Immediate. 

(4) Let Cy,?/] ClrW',u", v"]. Then, aziwi ~i a and q' = g" and 

u!w r z r p ~r u"wiv"w r z r p. Which is possible only if a ~ P ~ A and 
zi — z r — wi — w r — v" = □ and r' = u". But this is the case only 
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if C — (\,q,u\). And analogously for C[q f ,v!\ CrlW', u", ?/'], it 

may hold only if C — (A, g, uX). 

(5) The case (^ 1 ,^ 2 ) = (d^cK,) is immediate. Thus, let us assume 
ClrW, u r , v'\ Crr^'', u”, v"]. We see that this is possible only 

if a aziwiv"w r and q' = and u'wiv'w r z r (3 u"f3. Which 
holds only if a /? A and zi — z r = wi — w r — v' — v" = □ and 
u' = ?/'. That implies C = (A, g, rA). 

□ 

Lemma 54. Lei = 7(7(0 u, d x ),p, v, d 2 ). 

(1) There is a triple (q,WR,WR) G QxSxS, s.L ; for any (q',u f ) G QxS 
and d = we dat’e (7 = C q p ^f iff (p, p) = (g, Wd)- 

(2) For dx + d 2 € B*, C = cff/ 2 , iff C = (A, g, A) and u' = v = A 
and p — q. 


Proof. 


(i) 


( 2 ) 


Let us define Cr^ 7 ,?/] = (aziu',q f ,uwrZ2I3) (the configuration ob¬ 
tained by replacing q' and v! for the current state q and the re¬ 
verse right symbol wr of C) and Ci\q[,v!\ — ( az\WR,q f ,uu' Z 2 P) 
(the configuration obtained by replacing q' and v! for the current 
state q and the reverse right symbol wr of C), it is readily seen 
that Cyl v ff d — Cd\p,v] and, as a consequence, C — Cd\p,v] iff 
(P,v) = 


By direct computation, we see that Cp ^f R = (aziWRun'v,p, /3) and 


C 


jq' ,u' ,L 


_ qf 

p : v,R ~ {ot,p,vu r uwlzrP) . From which, we see that C — C pv f 

iff a = (3 = A and zr — zr — wr — wr — u — u r — v — U and 

„ -psq'u'.L . — . 

C = C pv R ltt Z L = Zr = W L = Wr = U = U = V = □. 


□ 


B.4. The adjoint of Um • We can now compute the adjoint of the operator 
Um • For this, we have already given the reverse transition 7 , but we also 
need to reverse the quantum transition function 5. For x G {0, «s, t}, let us 
take 

$X : % t 2 (s x x B x ) s.t. s x (p,v)(q,u,d) = 6 x (q,u)(p,v,d)* 
where Oq = B and = Bt = {i?}. Then, let us define 
£m = {(av,p,/3) e C M I (p,v) e %} 

—0 —s — t 

It is readily seen that and are a partition of £, since they 

—0 —s — t 

are pairwise disjoint and £ = U U Moreover, given 

C ~ (av R ,p,wv L /3) e Cm 
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we have that, (p,vr) G T x by definition, and when x = 0, that (p,vl) G 7o 
also (since by the definition of ground configuration, p G Qo U Qt implies 
that wvlP G E*). Thus, we can finally define 

W* M |C)= S x (P,Vd)(q,u,d) |C q , u ,d) 

(q,u,d)£S x xDj 

= $x(q,u){p,v d ,d)* |C q ud ) 

(g,lfc,d)&Sa;XDa; 

where C qwl = 7 S)Q (C, < 7 , it, rf). 

We remark that, w.r.t. the definition of in the range of the sum in 
we have S x x in the place of T x x B. In the case x = 0 , there is no dif¬ 
ference, since H) x = B and then we consider both the reverse deplacements; 
in the case x = s,t instead, B^ = {i?} and then we consider the R de¬ 
placement only. Technically, this is necessary as in these cases (p,vl) 0 7^, 
and then S(q,u)(p,VL, d) would not be defined. Indeed , for x = s,£, any 
configuration in £ X M can be entereded from an R deplacement only; then, in 
this case, to reverse the quantum transition functions it suffices to consider 
R reverse deplacements only. On the other hand, even in the definition of 
H 7Mi we might have restricted the sum to the R deplacement only, in the 
case x = s,t. Indeed, 

W M (|C»= S x (q,u)(p,v,d) |Cp,v,d) 

(p,v,d)eT x xB x 

for C — (cp g, u/3) G C^, and x = 0, s, t. 

Lemma 55. W ^ |C) G span(CB(£)) ; for every C G Cm- Then, defines 

an automorphism 

Wtf : span(CB(C)) -A span(CB(C)) 

of the linear space span(CB(C)). 

Proof. By case analysis, as in the proof of Proposition [5j □ 

We can also see that maps every <£ X M into C|^, which is indeed the 
converse of the fact that Wm maps £ X M into (Sf M . 

Lemma 56. 

(1) W M (span(CB(<^))) C span(CB(^)) 

( 2 ) WAf(sp an ( CB (^f))) C span(CB(£^)). 

Proof. By case analysis. □ 

We can now prove that defines the adjoint of Um- 

Lemma 57. The unique extension of to the Hilbert space ^ 2 (Cm) is 
the adjoint ofUM • 
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Proof. It suffices to prove that (W^ |D), |C)) = (D| Wm |C), for every C, D G 
Cm- Let C G £ x and D G with x,y G {0,s,£}. Since ( W^ |D) , |C)) = 
(D| Wm |C) = 0 if x 7 - y (by Lemma [56]), we have to analyse the cases x = y 
only. 

Let C ~ {aiWi.q.uPx) G (57, D ~ (a 2 VR,p,w 2 v L f3 2 ) G = 

7 ( 1 ?, g, r, d), and C PiVi d = 7 ((7,p, u, d). Let us start with the case x = 0. 


W |D), |C)) = 2 8 x (p,v d )(q',u',dy (D q / )U / )d I C) 

= ^ ^ ^ic((7 5 u )(Pj Vd-, d) (Dq'ju'jd 1^) 

= ^ ^ ^x((b u )(Pi Vdi d) (Dq jU? d | C) 

= XI s x(q,u)(p,v d ,d ) (D I Cp ;Vd;d ) 

= X dx(q,u)(p',v',d) (D | Cp/y id ) 

(p f ,v f ,d)eT x xB x 

= (D| W M |C) 

Where we have used the facts that (D q /, u /,d | C) = 0 if (q',u') 7 ^ (g,R), 
that (D I C p /y 5 d) = 0 if (p',v f ) 7 ^ (p,Vd) (by inspection, we see that in these 
cases the above configurations differ for the current state or for the current 
symbol), and that (D q?u? d | C) = (D | C p?Vd ,d ) 5 by Lemma 52 □ 


B.5. Unitarity of the time evolution operator. In the following, we 
shall complete the proof that the time evolution operator of a QTM is uni¬ 
tary. Firstly, we show that the time evolution operator Um of a pre-QTM 
M is an isometry iff the local unitary conditions hold (Lemma 58), then we 
shall see that, in the particular case of pre-QTMs, Um is unitary when it 
is is an isometry (Lemma 59). Therefore, a pre-QTM is a QTM iff its time 
evolution operator is unitary. 

Lemma 58. U^Um = 1 iff the local unitary conditions holds, that is, Um 
is an isometry iff M is a QTM. 

Proof. It suffices to prove that, for C G with x G {0, s, t}, 

W* m W m \C) = X S x (q,u)(p,v,d)W* M \ C p>Vid ) = 1 

( p,v,d)eTx xD 

iff the local unitary conditions hold. 

Let us assume 

C ~ (avRW,q,uzvLfd) G £ x 

C p ,v,d = l(C,p,v,d) c l'v/ =7(7 (C,p,v,d),q',u',d') 
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we have 


C p ,v,L ^ {av R , q, wvzvl/3 ) C PtVjR ~ ( av R wv , q, zv L j3) 


The cases x = s and x = t are immediate. For instance, if x — £, we have 
(<?,«) G Qt x E, and 5t(q,u)(p,v, d) = 1, when (p,v,d) — (g, R,i?), while 
5 t (q,u)(p,v,d) — 0, otherwise. Thus 


w* m Wm |C) = 


^M|C q ,u,R) = 


( q',u')eS t 


,u')(q,u,R)* 


f-qV-R 

'“q,Lj,R 


but 8 t (q',u')(p,u,R ) = 1, (p,u) 
wise; thus 

W* m Wm |C) 


= and 5 t (q',u')(p,u,r ) 



0, other- 


(by Lemma [53]) and analogously for the case x = s. Therefore, for x E {s, t}, 
since the equivalence holds for every pre-QTM. 

Let us now analyse the case x = 0. By definition, 


WL |C 


M |^-p,v,d/ 


= E = S o(q',u')(p,v' d ,,d') 

(iq',u',d') 


(-q'V,d' 

'“P^d 


for some v' d , that depend on C }pv j. By reorganizing the sums according to 
the cases d! = d and d' / d, we get 


P> v >d 


W* m Wm |C) = E S 0 (q,u)(p,v,d)W* M \C P:Vtd ) 

(p,v,d)eT x xB 

= E E S o(q,u)(p,v,d) S 0 (q',u')(p,v,d) 

(p,v,d)£Tx xB ( q',u')eS x 

+ E E S 0 (q,u)(p,v,L)5o(q',u')(p,v R ,R)* C^ R 

(p,v)eT x ( q'y)eS x 

+ E E s o(q,u)(p,v,R)5 0 (q',v!)(p,VL,L)* Cp V u R L 

{p,v)eTx {q',u')eS x 

By Lemma [53j this corresponds to 


W* m W m |C) 

= E \ s o{q,u){p,v,d)\ 2 |C> 

(p,v,d)eToxlb 

+ E E «)(p, v, d) d 0 (g', v , d )* |C[q', u']) 

(g / ,n / )e<So\{(g,w)} (p,^,d)e7bxD 

+ E E E do(q, u)(p, V, R) S 0 (q', u')(p, v L , L)* |C RL [q', u', v]> 

(q'y)eSo vgs peQouQt 

+ E E E u)(p,v, L)8 0 (q',u')(p,v R , R)* |C LR [q', u',v]) 

(q'y)eSo peQoUQt 
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where, C, C[q'u'], Clr'W,u' , p], and CrlW, n!, v] are never equal, does not 
depend on (p, v), and C, C[q'u'] are independent from d also. 

From which we may conclude that, W^Wm |C) = |C) for every C G (C°, 
iff for every q, g' G Qo U and p, u, v!vp , tr G £ 

1 = XI l^oU'OU'M)! 2 

(p,v,d)e7bxID> 

°= X s o(q,u)(p,v,d)5 0 (q , ,u')(p,v,d)* 

(p,v,d)eToX~D 

°= X s o(<h u )(PiV,R)$o(q',u')(p,v L ,L)* 

peQoUQt 

pCQoUQt 

that is, iff the local unitary conditions hold. □ 

Lemma 59. If Um is an isometry, then UmU^ = 1. As a consequence, 
Um is unitary. 

Proof. Since Um is an isometry, by Lemma [5lj we know that UmU^ is an or¬ 
thonormal projection and UmU* m — 1 iff (QUmU^ |C) = (C| WmWm IQ = 
1 for every C G (Cm- 
Let us assume that 

C ~ (av R ,p,wv L /3) e <t x 

C q ,u,d = 7(C q, u, d) Cp\vtd' = 7(7(0 q, U, d),p v' , d!) 

we have 

Cq,u,L - {mv r w, q, Up) Cq^R ~ (a, g, UWV L p) 


(C\W m W* m \C) 

= X s x(q,u){p,v d ,d)* (C| Wm |Cq iUi d) 

(q,u,d)£S x xD x 

= X s x(q,u)(p,v d ,d)*5 x (q,u)(p , ,v , ,d') 

(q,u,d)eS x xDcc 
(p' y ,d f )eT x x^x 


C I c 


^q,u,d 


p x , v ' ,d / 


By Lemma [54| we have that 
(1) when d — d! 


7^q,M 

'^p / ,v / ,d / 


r I 7^q,u,d 

L I L P 'y,d 


if (p', v') = (p, v d ) 
otherwise 
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(2) when d ^ d' 

/ _q,u,d \ f 1 ii C = (X,p, A) and p = p' 

[ 0 otherwise 

When x G {«§,£}, we have = {R} and then, by item 0. 

(C\W M W* M \C) = X \UQ,v)(p,v R ,R )\ 2 

(■ q,u)eS x 

Morevover, in both cases S s (q,u)(p,VR, R) = 0, with the exception of the 
case (q,u) = (p, vr), where 

(1) for x = 5, we have that vr = a G X and vr = a; 

(2) for x = £, we have that ^ = a G S and vr = a. 

Since 5 s (p, 5#)(p, rr, R) = 1, in both cases we get (C| Wm^m IQ = 1. 

Let us now consider the case x = 0. 

(C| WmWm IQ 

= XI l^o(g,w)(p,^,d )| 2 

(q,u,d)£So xDo 

+ X <*0(9, □ )(?, □, d)* S x (q, n)(p, □, d) (C | Cj’°’g) 

(q,d)£So xB 

where d — R when d — L. and d — L when d = R. We remark that the last 

addend may not be equal to 0 only when C = (A ,p, A) (see item <§ above). 

In order to analyse it, let us not consider a single configuration only, but 
a whole family of configurations 

C\p',v' R , v' L ] ~ (av' R ,p',wv' L f3) € £° 

which differ from C — (avR,p,wvR/3) G £° for the current state p and 
the current i?-reverse and L-reverse symbols vr and vr, respectively. More 
precisely, we take 

@C — {C[p\ v' R , v f L ] G £ | (p\v' R ,v' L ) gQx X 2 } 
which, by the definition of £° corresponds also to 

@C = {C[p\ Vr, Vr] I {p\v'r,v' l ) G (Qo LI Qt) x X 2 } 

From which, it is readily seen that 

\Bc\ = (|Qo U Qt\) |S | 2 
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Finally, let us take 

E (C'\WmW* m \C') 

C'eBc 

E (c[p'v'r,vl]|w |c[p'y R y L ]> 

O' ,v' R ,v' L )e{QoUQt)x'E 2 

= E E E \ 5 o(q,u)(p',v' R ,R)\ 2 

v' L eT, ( q,u)eS 0 (' p',v' R )eTo 

+ E E E M<huwy L ,L)\ 2 

v ' R eu ( q,u)eS 0 (p',v' L )eT 0 

+ E E <5 0 (g, □)(?/, □, 0* <$ 0 (g, □)(?>', □ , <0 (c | Cj' 

p'efiouQt ' 

(o , ,d)e(QoUQ s ) xD 

By the local unitary conditions, 

E <5o (q, □ )(?', □, 0* <So(9, □ )(?', □, d) = 0 

p'EQoUQt 

for every g G QoUQ s . Therefore, by taking into account that 5o(q, u)(p' , i?) 

does not depend on v f L , and that 5o(q,u)(p',v' L , L) does not depend on v ' R , 
we have 


E ( C '| |C'> = |£| E E l<^0 (q,u)(p',v' R ,R )\ 2 

CeB c ( q,u)eS 0 ( p',v' R )eTo 

+ |S| E E !*>(</, «)(P>LUI 2 

(g,u)e<S 0 (p',v' L )eTo 

= l s l E Y \ s o(q,u)(p',v',d )\ 2 

(q,u)eSo (p',v' ,d)e7bxD 


But by the local unitary conditions, Y2(p f v ' d)eToxB \S 0 (q,u){p',v',d)\ 2 = 1 
for every (q,u) G Sq. Thus, 

E <c'| WmWm |c') = |E| E 1 = I So U Qt\ |E| 2 = \B C \ 

CeB c ( p",v")eSo 
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Finally, let us recall that, by Lemma 


C G Cm- Therefore, for every C G C° 


51 


'M ’ 


(< C\W M W^\C ) < 1 for every 


1 > (C| W m W* m |C) 

= Y {t'\W M W* M \Q.') - Y (C'\W m W* m \C) 

C'eB c C'eB c \{C} 

= \B C \- Y < C 1 w mW* m |C') 

C'eB c \{C} 

>\Bc\-\B c \{C}\ = l 

That is, (C\W M W^\C) = 1. □ 

Theorem 60. A pre-QTM is a QTM iff its time evolution operator is uni¬ 
tary. 


Proof. If the time evolution operator Um of the pre-QTM M is unitary, and 
therefore an isometry, then the local unitary conditions hold (by Lemma [58]) , 
and thus M is a QTM. On the other hand, if M is a QTM, and therefore 
the unitary conditions hold, then Um is an isometry (by Lemma [58]) , and it 
is unitary indeed (by Lemma 59). □ 
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